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ABSTRACT 

Interstellar dust presents a significant challenge to extending parallax- 
determined distances of optically observed pulsational variables to larger vol- 
umes. Distance ladder work at mid-infrared wavebands, where dust effects are 
negligible and metallicity correlations are minimized, have been largely focused 
on few-epoch Cepheid studies. Here we present the first determination of mid- 
infrared period-luminosity (PL) relations of RR Lyrae stars from phase-resolved 
imaging using the preliminary data release of the Wide-Field Infrared Survey 
Explorer (WISE). We present a novel statistical framework to predict posterior 
distances of 76 well-observed RR Lyrae that uses the optically constructed prior 
distance moduli while simultaneously imposing a power-law PL relation to WISE- 
determined mean magnitudes. We find that the absolute magnitude in the bluest 
WISE filter is Mwi = (-0.421 ± 0.014) - (1.681 ± 0.147) logio(P/0.50118 day), 
with no evidence for a correlation with metallicity. Combining the results from 
the three bluest WISE filters, we find that a typical star in our sample has 
a distance measurement uncertainty of 0.97% (statistical) plus 1.17% (system- 
atic). We do not fundamentalize the periods of RRc stars to improve their fit 
to the relations. Taking the Hipparcos-derived mean ^-band magnitudes, we use 
the distance posteriors to determine a new optical metallicity-luminosity relation 
which we present in §5. The results of this analysis will soon be tested by HST 
parallax measurements and, eventually, with the Gaia astrometric mission. 

Subject headings: infrared: stars - methods: statistical - RR Lyrae: distance 
scale 
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1. Introduction 

RR Lyrae (RRL) pulsating variable stars are standardizable distance indicators at op- 
tical and near-infrared wavebands. In ^-band, their brightnesses are nearly standard, with 
a small metallicity dependence and deviation about {My) of ~0. 12-0. 15 mag (Hawlcy et al. 
1986; Fernley et al. 1998; Chaboyer 1999; Sandage & Tammann 2006). At near-infrared 
wavebands RRL brightnesses are a well-fit function of pulsation period, with an apparently 
negligible metallicity dependence (at K-band) and mean scatter from a period-luminosity 
(PL) relation of ~0.15 mag (Longmore et al. 1986; Sollima et al. 2006). The ability to in- 
fer distance to an RRL is chiefly restricted by the confidence in these empirically derived 
luminosity-metallicity and PL relations. 

There is good observational and theoretical motivation to believe that infrared photom- 
etry offers the ability to derive more tightly constrained PL relations for pulsational variable 
stars in general. It has been argued (Madore & Freedman 1998) and demonstrated (Frccd- 
man et al. 2008; Feast et al. 2008) that the scatter in these empirical relations is decreased 
at infrared wavelengths. Madore & Freedman (1998) cite the advantages: (1) The sensitivity 
of surface brightness to temperature is a steeply declining function of wavelength; (2) The 
interstellar extinction curve decreases as a function of increasing wavelength (being almost 
linear with 1/A at optical and near-infrared wavelengths); (3) At the temperatures typical 
of horizontal-branch stars, metallicity effects predominate in the UV, blue, and visual parts 
of the spectrum, where most of the line transitions occur, with declining effects at longer 
wavelengths. The overall insensitivity of infrared magnitudes of RRL, Cepheid, and Mira 
variables to each of these effects results in decreased amplitudes for individual pulsating 
variables, as well as a decreased scatter in the apparent PL relations. 

In this paper we present the first published mid-infrared PL relations for RRL variables. 
This is the first such work primarily because the requisite observational data has not existed 
previously. Since the farther reach of (brighter) Cepheid PL relations makes their study 
potentially more influential, the Spitzer Space Telescope has been used to derive mid-infrared 
PL relations for Galactic (Marengo et al. 2010) and Magellanic Clouds (Madore et al. 2009) 
Cepheids (the latter making use of SAGE survey data; Meixner et al. 2006; Madore et al. 
2009). These studies of Galactic (Magellanic Clouds) Cepheid mid- infrared PL relations 
reported best-fit dispersions of ~0.2 mag (~0.12 mag). 

RRL variables are particularly important local distance indicators because they are 
more numerous than Cepheids, and are observable within the Galactic disk and halo, within 
Galactic and some extragalactic globular clusters, and in the halos of neighboring dwarf 
galaxies (most notably, the LMC). Importantly, RRL variables can be used as stellar density 
tracers (e.g.. Vivas et al. 2001; Sesar et al. 2010) to map the structure of stellar distributions. 
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In this article we derive mid-infrared RRL PL relations by analyzing observations of 
76 RRL-type stars conducted with the Wide-field Infrared Survey Explorer (WISE) satellite 
(Wright et al. 2010) and made available through the preliminary data release of the first 105 
days of science data^. We use a modified Lomb-Scargle (Lomb 1976; Earning 1963; Scarglc 
1982) period-finding algorithm to calculate the pulsation periods from both the WISE data 
and the very well-observed Hipparcos light curves of the same sources. We derive mean 
flux-weighted WISE magnitudes from the best-fit harmonic models of this Lomb-Scargle 
analysis; these observed magnitudes, along with the Hipparcos estimated periods, are used 
to estimate the WISE PL relations. The actual PL fitting is conducted through a Bayesian 
approach using a priori distance information and simultaneously fits the Wl, W2, and W3 
PL relations. Our methods have general applicability, and can be used to robustly fit PL 
relations at other spectral wavebands. Our resulting mid-infrared PL relations are tightly 
constrained with absolute magnitude prediction uncertainties as small as 0.016, 0.016, and 
0.076 mag at 3.4, 4.6, and 12 /zm, respectively. 

The paper is outlined as follows. We present a brief description of the WISE and 
ancillary data in §2, followed by an explanation of our analysis methods in §3. (In §A we 
demonstrate the viability of period recovery with WISE data, and highlight the potential 
for discovery of new RRL variables and other short-period variables with the WISE single 
exposure database.) We describe the Bayesian method of deriving mid- infrared PL relations 
in §4 and discuss the results in §5. Finally, we present conclusions in §6. 

2. Data Description 

WISE has imaging capabilities in four mid-infrared bands: Wl centered at 3.4 /xm, W2 
at 4.6 /xm, W3 at 12 fim, and W4 at 22 /im. The satellite is in a polar orbit and scans the 
sky in great circles with a center located at the Sun and with a precession rate of 180° every 
six months (Wright et al. 2010). WISE completes about 15 orbits a day and the field of 
view of the detectors is 47 arcmin on a side. This configuration allowed WISE to scan the 
entire sky in six months, with a minimum of 8 (median 12) single-frame exposures. Sources 
near ecliptic poles receive the most repeat coverage in time and sources near the ecliptic 
have the smallest number of observed epochs. WISE was launched on 2009 December 14 
and operated until 2011 February 17. This mission duration provided two full scans of the 
sky. However, the hydrogen coolant ran out in 2010 October, halting data acquisition in the 
W3 and W4 bands. 



http : //wise2 . ipac . caltech.edu/docs/release/prelim/ 
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The present study was conducted using data from the single-exposure database of 
the WISE prehminary data release, which was made public on 2011 April 14 through the 
NASA/IPAC Infrared Science Archive^. The preliminary data release includes the first 105 
days of mission data and covers about 57% of the sky. 

The catalog of RRL variables used in the present study is derived from work by Fernley 
et al. (1998). The catalog contains 144 relatively local (<2.5 kpc) RRL variables selected 
from the Hipparcos catalog (M.A.C. Ferryman & ESA 1997) with color excess and metallicity 
measurements from previous literature. Of the 144 RRL variables in our starting catalog, 
77 were associated with sources in the WISE preliminary data release. We reject one light 
curve (V*EZLyr) because the reported WISE photometry does not indicate a periodic source 
(using the Hipparcos period; the source is also a strong outlier from our FL relation fits). 
All target sources, save the prototype RRL itself (V*RRLyr), are too faint for WISE to 
produce reliable W4 photometry and so we must ignore the longest wavelength data in the 
subsequent analysis. 

To perform a proper analysis of the mid-infrared FL relation, in addition to the periods 
and observed WISE magnitudes of each RRL, we also need a prior guess of the distance to 
each object. Here, we describe how we determine a prior distribution on the distance modulus 
of each RRL. First, we compute the Hipparcos periods and mean flux- weighted magnitudes 
(mhip) using the same Lomb-Scargle based methods that we apply to the WISE light curve 
data (see §3). Discrepancies with the published periods of Fernley ct al. (1998) were minimal 
(see Appendix A). Unlike Sollima ct al. (2006), we do not find that the periods of RRc type 
RRL variables need to be fundamentalized by adding a constant term (Alog^g-^ ~ 0.13) in 
order to improve the FL relation scatter. Following Gould & Fopowski (1998) we determine 
values of the apparent Johnson V^-band magnitude (my) and the effective extinction, Ay^^s, 
that differ slightly from Fernley et al. (1998). This is achieved by making use of the line-of- 
sight extinction from the Schlegel et al. (1998) (SFD) dust maps and by assuming a Galactic 
scale-height model for the dust (such that not all SFD dust lies in between us and the RRL). 
In particular, we determine an effective extinction for the ith RRL in the sample as: 

E{B - r)eff,i = E{B - V)sFD,i (1 - eM-M/h]) , (1) 

where E{B — V^)sFD,i is the differential SFD extinction towards source i, Zi is the scale height 
of ith source above the Galactic plane'^, and h = 130 pc is the exponential scale height 



■^http : / / irsa . ipac . caltech . edu/ 

■^Note that Zi will depend on the value of distance determined, however, we simply use the coordinate 
information provided by Maintz & de Boer (2005) when available or transform the sky coordinates using the 
Fernley et al. (1998) distance results when necessary. Our results are not sensitive to the precise value of z^. 
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assumed for the dust in the disk (Gould & Popowski 1998). To convert mhip to my we adopt 
the prescription from Gould & Popowski (1998): 

mv,^,es = mhip.i -X- 0.2E{B - VU;, (2) 

where X = 0.09 for RRab types and X = 0.06 for RRc types. We assume a 15% error on 
E{B — V)e{i^i- Finally, we produce an extinction corrected magnitude uiy^ = mv,i,es — Av,es,i, 
where Av,es,i = 3.1 x E{B — V)es,i, and the factor R = Ay/E{B — V) = 3.1 from Schultz 
& Wiemer (1975). These values of extinction-corrected l^-band magnitudes (and associated 
errors) are reported in Table 1. 

To determine the prior distance modulus, /io,j, for the zth source, we need a prescription 
for determining the absolute l^-band magnitude given the metallicity of the star (there is 
no known relationship of period and luminosity at V^-band for RRL variables). We adopt 
the My-metallicity relation given in Chaboyer (1999), where we use the metallicity data as 
provided in Fernley et al. (1998). Explicitly, the My— [Fe/H] relation used is 

Mv = (0.23 ± 0.04)([Fe/H] + 1.6) + (0.56 ± 0.12). (3) 

The calculated values of My j for each source are given in Table 1. 

Finally, we compute the prior mean of the distance modulus of the ith RRL as fio,i = 
rriy- — My^i, with the uncertainty in this quantity propagated assuming the errors on ■ 
and Myj are Gaussian and uncorrected. The values of /io,i (and cr^g J (Table 1) represent 
our best estimates of the distances (and errors) using the body of work on RRL variables 
at visual bands prior to analyzing the WISE data and prior to the Hubble Space Telescope 
(HST) parallax result on V*RRLyrae itself. Note however that our prior estimate on distance 
modulus to V*RRLyrae (/xq = 7.042±0.125) is consistent with that found directly from HST 
parallax measurements (/i = 7.090 ± 0.063; Benedict et al. 2002). 



3. Light Curve Analysis Methods 

In order to determine the PL relations for WISE, we need to calculate an apparent 
brightness. Following common practice (Fernley et al. 1998; Liu & Janes 1990), we define 
the brightness of each source as the mean flux, converted to a magnitude. As we expect 
possible poor phase sampling in the WISE data, we use a model — based on a modified 
Lomb-Scargle algorithm (Richards et al. 2011), which allows for data uncertainty and a mean 
flux offset — instead of the observed data points, to determine this mean (see Figure 1 for an 
example). At signiflcant peaks in the periodogram, this model construction attempts to flt 
as many as 8 harmonic components — at frequencies which are multiples of the fundamental 




Fig. 1. — Hipparcos and WISE light curves of V*UPic with fitted models. The solid horizon- 
tal line is at the mean flux magnitude and the dashed lines represent the uncertainty of this 
mean magnitude. The model fitted to the Hipparcos (WISE Wl) data uses 8 (7) harmonics. 
V*UPic is our best observed WISE source with 35 observations over ~2.3 days. The data is 
phase-shifted so that the minimum of the model in each band is at phase = 1. 
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frequency — in addition to the fundamental frequency component. Complex models are 
penalized using generalized cross validation (e.g., Hastie et al. 2009) to prevent over-fitting. 
The resulting model curves are smooth, typically dominated by the presence of 4-6 harmonics 
for Hipparcos data, and can be used to calculate the flux integral and its uncertainty. For 
the case of Hipparcos, we find that the difference between our flux estimates and those from 
Fernley et al. (1998) exhibit an rms scatter of 1.4% with no systematic difference. 

Applied directly to the WISE data, our period finding framework accurately recovers the 
majority of the RRL periods directly from the WISE data (see Appendix A). For all WISE 
mean-magnitude estimates, we force the Lomb-Scargle model to use the best-fit Hipparcos 
periods as the fundamental frequency. We note that our mean-magnitude estimates remain 
unchanged, within their uncertainties, if we instead use the best-fit WISE periods. 

4. Deriving the PL Relations 

Using the full sample of 76 WISE RRL variables, we derive the empirical PL relationship 
in each of the bands Wl, W2, and W3. For each RRL in the sample, we estimate the 
observed magnitude, niij, and pulsational period. Pi, using the methods outlined in §3. 
Here, i = 1, ...,n indexes the RRL variables and j = 1, 2, 3 indexes the WISE bands. 

Our statistical model of the PL relationship is^ 

rriij = fJ'i + Mo J + aj logio(Pi/Po) + Uj, (4) 

where /ij is the distance modulus for ith RRL, Mqj is the absolute magnitude zero point for 
the jth WISE band at P = Pq, where Pq = 0.50118 day is the mean period of the sample, 
and aj is the slope of the PL relationship in the jth band. We assume that any extinction 
is negligible in these bands. The error terms etj are independent zero-mean Gaussian ran- 
dom deviates with variance (aam.^)^, which describe the intrinsic scatter in the rriij about 
the model, where a is a free parameter which is an unknown scale factor on the known 
measurement errors, Crm/^- We fit the model (eq. 4) using a Bayesian procedure, described 
below. 

A Bayesian approach to this problem is appropriate because for each RRL we have a 
priori distance information from previous V-band RRL studies. For each RRL in our sample. 



^In principle, there could be a metallicity dependence, but we found such a dependence was negligible in 
the WISE bands. See §5. 

^The average measurement error, tT,„, is 0.013, 0.013, and 0.045 mag in Wl, W2, and W3, respectively. 
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we determine a prior on its distance modulus using the steps outlined in §2. For the star 
V*RRLyr, we adopt the HST distance estimate of Benedict et al. (2002) as our prior. These 
priors encompass the full amount of information that we have about each source's distance 
before looking at the WISE data. The key in our analysis is that while the distance to any 
RRL could be changed within its prior to fit a perfect PL relation in a single band, the 
simultaneous fitting of a power-law PL relation in all bands (with as little intrinsic scatter 
as possible) tightly constrains the distance of each source. Bayesian fitting of the PL model 
allows us to obtain: 

• posterior distributions on the distance to each RRL, given the WISE data, 

• posterior distributions on the absolute magnitude zero point and slope of the PL rela- 
tionship in each WISE band, and 

• an estimate of the amount of intrinsic spread of the data around the PL relationship. 

The end goal, of course, is to use the estimated PL relationship to accurately predict the 
distance to each newly observed RRL from its period and observed WISE light curve. Fur- 
thermore, we want to make these predictions with an accurate notion of the amount of error 
in each predicted distance, as those errors will propagate to subsequent studies. 

Bayesian fitting of linear models is thoroughly described in Gclman ct al. (2003). Here, 
we summarize our procedure for analysis of the WISE PL relationship. First, we assume 
a normal (Gaussian) prior distribution on each of the distance moduli with mean /io,i and 
standard deviation cr^g^, as described above. For the other parameters in our model (eq. 4), 
we assume a flat, noninformative prior distribution. For convenience, we rewrite the model in 
matrix form as m = X/3 + e, where m is a vector of the 3n measured WISE mean-magnitudes, 
/3 is a vector of the n + 6 parameters (/i,Mo,a) in the PL model, X is the appropriate 3n 
by n + 6 design matrix for eq. 4, and e is a vector of the zero-mean, normally distributed 
random errors with covariance matrix cr^diag(cr^). 

Including an informative prior on fi is equivalent to adding extra prior "data points" to 
the analysis. In our model, these "data points" are /io,j = /^i + cT/xoi^i) where is a normal 
random variate with mean and variance 1. This prior information on fi induces the model 
= X*/5 + e*, where 




(5) 




(6) 
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aMiag(a^) 



0„,3„ diag((j 



2 
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~ A^(0, S*), and denotes the multivariate normal distribution. Here, indicates the 
n X n identity matrix and Om,n is the m x n matrix of Os. 

Posterior distributions for the parameters of interest can be derived in a straightforward 
manner using the entities in eqs. 5-7. The joint posterior distribution, P{/3,a\m, P), can 
be sampled by first drawing from P{a\m,P) and then, conditional on that draw, selecting 
from P{P\m, P,a). The posterior distribution for (3, conditional on the value of a, follows 
the multivariate normal distribution, 

/3|m,P,a~iV(^,(X:S;iX,)-i) (8) 

where /3 is the standard maximum likelihood (weighted least squares) solution, 

5=(x:s;X)-Xs;V. (9) 

Unlike the posterior distribution of P (given a), the posterior distribution of a does not 
follow a simple conjugate distribution. Instead, the distribution follows the form 

where the prior on /3 is proportional to the informative prior on /i, the flat prior on a is 
P(o"^) oc a~^, and the data likelihood L is the product, over all observed magnitudes, of the 
Gaussian likelihood of the data given the model (eq. 4) with all parameters specified. 

We draw samples from our joint posterior distribution P{/3, a\m, P) using eqs. 8 and 10 
in conjunction. In practice, we compute*^ P{a'^\m, P) over a fine grid of a values using eq. 10, 
and then draw a sample of a from this density. For each sampled a, we subsequently draw 
a (3 from eq. 8, conditional on the drawn a value. We repeat this process 10,000 times to 
characterize the joint posterior distribution. Using a large sample from this joint posterior 
distribution, we can compute quantities of interest such as the maximum a posteriori slopes 
and zero points of the PL relationship of each WISE band, the intrinsic scatter of the data 
around the PL relationship in each band, and the spread in the a posteriori distribution of 
the PL parameters (see Figs. 2-4). 

^Assuming that (3 = f3. Several iterations show that the posterior distribution of a is insensitive to the 
assumed choice of /3. 
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5. 



PL Relations Discussion 



Bayesian analysis of the WISE RRL variables shows a strong PL relationship in each 
of the three bands. The maximum a posteriori estimates (and corresponding errors) of the 
slopes and absolute magnitude zero points for each of the three bands (eq. 4) and the joint 
posterior distributions of these parameters are plotted in Figures 2, 3, and 4 (see also §6). 
At the mean period (Pq = 0.50118 day) of the sample, we achieve an absolute magnitude 
prediction error of 0.016, 0.016, and 0.076 mag in Wl, W2, and W3, respectively. Therefore, 
for an RRL of period near 0.5 day observed in WISE Wl or W2 bands, we can predict the 
absolute magnitude of that object to within 0.016 mag, which corresponds to a fractional 
distance error of 0.7%. 

The width of the absolute magnitude prediction bands becomes slightly larger as one 
moves to larger or smaller periods, as the model is less constrained in those regions. How- 
ever, the prediction uncertainty remains low throughout the full period range, even at the 
extremes. For example, at a period of 0.3 day, the absolute magnitude prediction error is 
0.037, 0.037, and 0.084 mag (1.7, 1.7, and 3.9% fractional distance error); at a period of 0.7 
day, the prediction errors are 0.026, 0.026, and 0.079 mag (1.2, 1.2, and 3.6% fractional dis- 
tance error) for Wl, W2, and W3, respectively. Figure 5 plots the estimated PL relationship 
in each WISE band, plus the ilex prediction intervals. For each newly observed RRL, the 
true absolute magnitude is expected to reside within the prediction interval. 

Along with estimating the PL relationship for each band, our fitting procedure supplies 
a posterior distribution for the distances of each of the RRL in our sample. In Table 2 
we report the posterior means along with the 68% and 95% posterior credible sets for the 
distance to each of the 76 RRL used to fit the PL relationships. We also list the separation 
between prior and posterior distance moduli in units of a, defined as 



where /2 and o"^ denote the means and standard deviations of the distributions, respectively. 
We note that there is, for most RRL variables, a close correspondence between the prior and 
posterior distances, as |A| < 2 for all but 2 sources in our sample (V*ANSer and V*IIKPup). 
Figure 6 shows a plot of prior versus posterior distance moduli, including a residuals plot, 
which shows again that, within their errors, the posterior distance distributions are consistent 
with the prior distance distributions for almost all the RRL. 

Recall that for V*RRLyr we use the well-measured HST parallax result, which corre- 
sponds to 262 ± 7.5 pc. Our posterior fit distance for V*RRLyr is 253 ± 2 pc, which is 
consistent with the HST distance at a level of 1.2o". We also get a consistent prediction for 
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Slope (a) Count 

Fig. 2. — Contour plot and histograms of 10,000 samples from the posterior distribution of 
the slope (ai) and absolute magnitude zero point (Mo,i) of the period-luminosity relation 
for Wl. Our data constrain ai to -1.681 ±0.147 and Mo,i to -0.421 ±0.014, with negligible 
correlation between those parameters. Levels in the 2D contour plot are at the 99.9, 99, 
97.5, 95, 90, 85, 80, and 70th percentile. 
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Fig. 3. — Same as in Figure 2, for W2. Our data constrain 02 to —1.715 ± 0.147 and Mo,2 
to —0.423 ± 0.014, with negligible correlation between those parameters. 
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Fig. 4. — Same as in Figure 2, for W3. Our data constrain as to —1.688 ± 0.150 and Mo,3 
to —0.493 ± 0.015, with negligible correlation between those parameters. 
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Fig. 5. — Period-luminosity relations for Wl, W2, and W3 (left to right), as found by 
our Bayesian fitting method. In each figure, the solid line shows our model's prediction 
of the RRL absolute magnitude, as a function of RRL period. The dashed lines show the 
±lcr prediction intervals; we expect the true absolute magnitude to fall within the dashed 
fines for 68% of all newly observed RRL variables. The top panel of each plot shows the 
residual spread around the best fit model, showing the small variance, cr^, in the intrinsic 
scatter around the PLRs. In the figure, RRab are plotted as circles and RRc as squares. 
Blazhko-affected RRL are indicated by filled points. 
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Prior ij 

Fig. 6. — Prior versus posterior distance moduli (bottom) and residual difference (top). 
The posterior distance modulus distribution is determined by sampling from the density in 
eq. 8, which considers the evidence in all three WISE bands as well as the prior distance 
information. As is evident, the posterior distances are consistent with the prior distance 
distribution, within their errors. In the figure, RRab are plotted as circles and RRc as 
squares. Blazhko-affected RRL are indicated by filled points. 
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V*RRLyr if we do not use V*RRLyr itself in the PL analysis (Fig. 7). That our analysis for 
the source with the most highly constrained distance prior is consistent with those results is 
further evidence of its accuracy and apphcabihty. 

To check the sensitivity of our results to the prior distances used, we analyze the changes 
in our posterior distance estimates under systematic prior offsets. We first note that the 
prior n estimate for V*RRLyr using the HST parallax result differs by 0.05 dex from the 
Hipparcos V-band estimate. To estimate the amount of systematic error in our posterior 
distance estimates, we inflate the prior mean distance modulus by 0.05 dex for a random 
50% of the RRL before running our Bayesian PL model fitting. As a result, the posterior 
distance moduli increase by an average of 0.023 dex, implying a systematic error of 1.17% 
on distance estimation. We take this to be a reasonable estimate of the systematic error. 

As a further sanity check, in Figure 7 we compare the prior, posterior, and prediction 
/i densities for a few RRL variables. The prediction density for each RRL was computed 
by holding out that particular RRL during the model fitting, and then applying the fitted 
model to predict the distance modulus of that source. We find that these "cross-vahdated" 
fi prediction densities are very consistent with the posterior densities, suggesting that 
the model is stable and that small changes in the set of RRL used to fit the model do not 
cause any substantial differences in the model. Furthermore, those densities are much more 
narrow than the prior densities, showing that the WISE data can constrain the distances 
to a great degree. Additionally, we see that both the posterior and prediction densities fall 
within high-probability regions of the prior distribution for three of the four stars, meaning 
that our model is in good agreement with the prior distances. Note that the one discrepant 
star plotted, V*ANSer, has the second largest discrepancy between prior and posterior /j, 
densities, after V*HKPup (Table 2). 

We also test whether including RRL metallicity into the model improves the PL rela- 
tionship fits. To do this, we add an additional term, jjZi, to our model (4), where Zi is the 
metallicity of RRL i and jj is the slope of the magnitude-metallicity relationship for the jth 
WISE band. Fitting this model, we find that jj has a significantly positive value, but that 
the predictive power of the new model, as measured by the width of the prediction intervals 
around the absolute magnitudes, does not differ from the original model which neglected 
metallicity. Furthermore, if we first subtract from the absolute magnitudes the fit of the 
model that uses only period, we find no relationship between the residuals and metallicity 
(slope of —0.00034 ± 0.00151). Including only period in the model achieves significantly 
better fits than including only metallicity, with half as much residual scatter. These results 
suggest that all of the absolute magnitude information encoded in [Fe/H] is already contained 
in the period, and so metallicity need not be added as a covariate in the model. 



V* AEBoo 
P =0.31 days 




V* ANSer 
P =0.52 days 



V* RRLyr 
P =0.57 days 




V* SSOct 
P =0.62 days 




Fig. 7. — Comparison of the prior, posterior, and prediction density for the distance modulus, 
/i, for 4 RRL in our sample. From top to bottom, these sources are V*AEBoo, V*ANSer, 
V*RRLyr, and V*SSOct. In each plot, the broad black curve represents the prior fi density, 
the red (hidden) curve is the posterior density of fi, and the green (foreground) curve is the 
prediction density for /i, which is found by holding that source out during model fitting and 
then predicting its fi with the built model. For three of the four stars, the posterior densities 
are in good agreement with the prior; V* ANSer is the second most discrepant star in our 
sample (after V*HKPup). In all cases, the posterior and predictive densities are much more 
precise than the prior densities and are in very close agreement to one another. 
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Finally, we derive an empirical My — [Fe/H] relationship using the posterior mean fi 
values from our Bayesian fitting to the WISE data. From this data, our best fit relationship 
is Mv = (0.10 ± 0.02) ([Fe/H] + 1.6) + (0.59 ± 0.10), which differs significantly in its slope, 
but not in its intercept value, to the Cliaboycr (1999) relationship — in eq. 3 — that was used 
to compute the original distance priors. Figure 8 shows a scatterplot of our estimated My as 
a function of metallicity; there is significant scatter around the empirical relationship, with 
a handful of large outliers. We also overplot the Chaboyer (1999) relation to demonstrate 
that both relations fit the data reasonably well. We qualify our new My — [Fe/H] relation by 
noting that the relatively constrained metallicity range of our sample RRL variables limits 
the relation's applicability at other metallicities. As RRL metallicity deviates from ~— 1.5 
the uncertainty in the slope of the relation rises steeply. 

At a glance, the nontrivial difference between the Chaboyer (1999) My — [Fe/H] relation- 
ship used to calculate our distance priors and the new My— [Fe/H] relationship we derive 
using our distance posteriors could indicate an inconsistency in the Bayesian approach to 
our PL relation fits. In particular, this discrepancy may suggest that the large spread in the 
prior distance distributions has allowed the Bayesian fitting technique too much freedom in 
computing posterior distances. To test this, we run a simple weighted least squares regression 
to fit each of the PL relations, fixing the distances at the exact values from the Chaboyer 
(1999) My— [Fe/H] relationship (without using a Bayesian fitting method to update the dis- 
tance estimates). This simpler fitting method results in statistically identical slope and zero 
point parameters for all three WISE bands. The scatter about the least squares fit, however, 
increases to 0.12 mag in Wl and W2, and 0.15 mag in W3 (from 0.016, and 0.076 from the 
Bayesian method). This increased scatter is expected, since the primary purpose of applying 
the Bayesian fitting technique is to reduce this scatter by simultaneously finding more accu- 
rate distances through the posterior distribution (i.e., updating the distance estimates given 
the WISE data). We can thus state confidently that the discrepancy between the Chaboyer 
(1999) My — [Fe/H] relationship and the new My — [Fe/H] relationship that we derive does 
not affect the PL relation fits. 



We have presented the first calibration of the RRL period-luminosity relations at three 
mid-infrared wavelengths. Our estimated PL relations, tied to the Vega magnitude system, 
are: 



6. 



Conclusions 



Mwi 

Mw2 



(-0.421 ± 0.014) - (1.681 ± 0.147) logio(P/0.50118 day) 
(-0.423 ± 0.014) - (1.715 ± 0.147) logio(P/0.50118 day) 



(11) 
(12) 
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-2.5 -2 -1.5 -1 -0.5 

[Fe/H] 

Fig. 8. — My — [Fe/H] relation as derived from our new, WISE-constrained RRL dis- 
tances. In the figure, RRab are plotted as circles and RRc as squares. Blazhko-affected 
RRL are indicated by filled points. The black thick lines depict our fit of the relation. 
My = 0.10([Fe/H] + 1.6) + 0.59, while the amount of intrinsic scatter about the fit is ±0.10 
mag, represented by the black dashed lines. The blue thinner lines depict the relation from 
Chaboycr (1999), My = (0.23 ± 0.04)([Fe/H] + 1.6) + (0.56 ± 0.12), with dash-dot lines 
showing the ±la bounds of the relation. 
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Mw3 = (-0.493 ± 0.015) - (1.688 ± 0.150) logio(P/0.50118 day). (13) 

These relations achieve an absolute magnitude prediction error as low as ±0.016 mag in WISE 
bands Wl and W2 (rising to 0.076 mag in W3) near the mean period value Pq = 0.50118 
day. Using these relations we calculated new distances to our sample of RRL stars with a 
mean fractional distance error of 0.97% (statistical) and 1.17% (systematic). 

We further demonstrated that the posterior distances resulting from the newly-derived 
PL relations are consistent with the prior distance distributions. An attempt to find an 
independent, statistically significant metallicity dependence in the mid-infrared PL relations 
confirmed the mid-infrared relations' independence from metallicity effects. Additionally, we 
applied our posterior distance estimates of our 76 RRL sample to fit a new absolute \^-band 
luminosity-metallicity relation. 

Perhaps the most significant contribution possible of the RRL PL relation is a well- 
constrained measurement of the LMC distance. The distance modulus of the LMC is a hugely 
consequential value in the extension of the distance ladder out to cosmological scales, and 
the the subsequent calculation of the Hubble constant, Hq (Schaefer 2008). The mid-infrared 
PL relations presented here will allow future studies of LMC RRL variables conducted with 
Spitzer (warm) or possibly the James Webb Space Telescope (JWST) to measure reliable 
LMC distances with error at the ~2% level or lower^. It is conceivable that a comprehen- 
sive mid-infrared survey of LMC RRL variables would enable the three-dimensional stellar 
structure mapping of the LMC with ~1 kpc resolution. 

The accuracy of any estimate of the PL relation is influenced by the accuracy of the a 
priori distances for the RRL sample used. Soon the HST parallax measurements (Benedict 
2008) of V*RZCep, V*UVOct, V*SUDra and V*XZCyg will be pubhshed. Our results in 
Table 2 serve as predictions of what will be found for the first three of those sources (once the 
WISE data on V*XZCyg is released, eqs. 13 could be used to postdict the HST result). The 
Gaia satellite of the European Space Agency, a 5-year astrometry mission to be launched in 
mid-2013, promises trigonometric parallax measurements of all field RRL variables within 3 
kpc with individual accuracy cr(7r) < 3% (Cacciari 2009). Although these measurements will 
not be available for many years to come, they have tremendous potential to further constrain 
the PL relations presented herein. In doing so, we can hope to study Galactic substructure 
well into the optically-obscured Galactic plane and further improve the resulting distance 



^Note that the current absolute cahbration uncertainty of WISE relative to Spitzer is 2.4, 2.8, 4.5% 
(W1,W2,W3, respectively), as provided in the Explanatory Supplement to the WISE Preliminary Data Re- 
lease Products — http://wise2.ipac.caltech.edu/docs/release/prelim/expsup/sec4_3g.html. This 
would dominate over the errors in our WISE-determined distance measure. 
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estimates for the LMC and beyond. 
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A. WISE Period Recovery 

There are two primary concerns with using WISE data to discover short-period variable 
stars. First, the number of observations on any given patch of sky is small (minimum 16 
in the final WISE dataset) and determined primarily by the ecliptic latitude. Secondly, the 
peak-to-trough amplitude of pulsating variables is significantly decreased at mid-infrared 
wavelengths as compared to optical wavelengths (~0.2 mag in Wl compared to ~1 mag in 
V for RRL variables). Our analysis shows that even with these disadvantageous factors, 
the WISE light curves can yield accurate periods quite often. Peaks in the periodogram are 
expected to have frequency widths ~ 1/T, where T is the time spanned by the observations. 
We note that our best-fit frequencies, determined on a grid of frequency steps 0.01/T, agree 
with well with those of Fernley et al. (1998) (to better than 0.2/T typically). We plot an 
example periodogram using the W2 light curve of an RRL with the median number of WISE 
observations (14) in Figure 9. 

For the fitting of PL relations, it is important to have accurate log-Period estimates 
or, equivalently, accurate fractional period estimates. We describe the (in) accuracy of a 
recovered period by the simple fractional error as compared to the known, true period. 

Recovered Period Fract. Error = iPm — Pt| /Pt (Al) 

with Pin the period measured solely from the WISE light curve and Pt the true period as 
measured from the Hipparcos light curve. 
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Period [d] 



Fig. 9. — Periodogram generated from the 14 epochs of W2 data for V*MSAra. To envelope 
the reasonable period range of RRL we plot from 0.2 to 1.2 day. The archival Hipparcos 
period of 0.525 day is well recovered from the WISE data (peak at 0.522 day). 
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Figure 10 shows histograms of the period recovery accuracy for each WISE band relative 
to Hipparcos and illustrates that nearly all WISE light curves produce accurate periods in 
the shorter wavelength bands Wl and W2. 

To explore how the number of epochs in a light curve affects period recovery, we plot 
in Figure 11 recovered period fractional error as a function of the number of observations 
for band W2. Because of the survey strategy of WISE, a larger number of observations 
typically indicates both increased temporal resolution (increased frequency of observation) 
and increased total light curve timespan (duration between first and last observation). As 
expected, there is a general trend of reduced recovered period fractional error with increasing 
number of observations. Beyond 20 observations, the typical period error is < 2%. 

Any period-finding algorithm must distinguish a shape for the light curve. That is, a 
phased light curve must be smoothly varying in that the uncertainty in the brightness at any 
phase point is considerably smaller than the amplitude of the light curve. As the photometric 
uncertainty increases relative to the amplitude, there is an effect of "vertical smudging" in 
which the light curve shape becomes less distinguishable. The flux amplitudes of RRL 
variables are about two times smaller in the mid-infrared as compared to the visual band. 
To investigate if this plays a factor in period recovery with WISE light curve data. Figure 
12 plots recovered period fractional error as a function of light curve amplitude/ (mag error) 
for band W2. Although we would expect to observe decreased period error with increased 
amplitude/ (mag error), this is not observed. We can conclude that "vertical smudging" is 
at most a non-dominant source of error in the recovered periods. 
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Fig. 10. — Distribution of recovered period fractional error across the three WISE bands. 
The plots are cropped to a maximum fractional error of 0.5 and the resultant number of 
excluded light curves is noted in each plot. 
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Fig. 11. — Recovered period fractional error plotted as a function of number of WISE 
observations for band W2. As expected, there is a general trend of reduced recovered period 
fractional error with increasing number of observations. Beyond 20 observations, the typical 
period error is < 2%. 
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Fig. 12. — Recovered period fractional error plotted as a function of light curve 
amplitude/ (mag error) for band W2. Although we would expect reduced recovered period 
fractional error with increasing amplitude/ (mag error), this trend is not obvious in the plot. 
"Vertical light curve smudging" is not a significant source of error in the recovered period 
for our dataset. 
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Table 2. RRL Distance and My Posteriors from Bayesian Analysis 



Name 


dbest ^ 


[d - 1(T, d + la] 


[d - 2(T, d + 2a] 


^(Mprior Mpost) ^ 




My 


c 




[pc] 




[pc] 


[pc] 


[No. of a] 




[mag] 


V*AACMi 


1330.7 


[1320 


4, 


1341. 1[ 


[1310.1 


1351.6] 


1.27 





716±0 


020 


V*AEBoo 


914.9 


[900 


6, 


929. 5[ 


[886.6 


944.2] 


-1.25 





770±0 


035 


V*AFVir 


1643.8 


[1631 


2, 


1656. 6[ 


[1618.6 


1669.4] 


-0.20 





647±0 


021 


V*AMVir 


1360.2 


[1347 


3, 


1373.2] 


[1334.5 


1386.4] 


-0.09 





625±0 


023 


V*ANSer 


1120.8 


[1112 


3, 


1129.3] 


[1103.8 


1137.9] 


2.44 





569±0 


019 


V*APSer 


1154.9 


[1139 


2, 


1170.8] 


[1123.6 


1187.0] 


-0.81 





668±0 


031 


V*ARHer 


1344.8 


[1334 


4, 


1355.4] 


[1324.0 


1366.0] 


0.20 





604±0 


019 


V*ARPer 


597.1 


[591 


8, 


602.4] 


[586.6 


607.7] 


-1.24 


1 


029±0 


024 


V*ATScr 


1499.6 


[1479 


3, 


1520.3] 


[1459.2 


1541.2] 


-0.11 





476±0 


032 


V*AUVir 


1515.7 


[1495 


2, 


1536.5] 


[1475.0 


1557.6] 


-0.22 





611±0 


031 


V*BBEri 


1389.6 


[1378 


4, 


1400.9] 


[1367.3 


1412.4] 


-0.19 





648±0 


020 


V*BCDra 


1424.1 


[1406 


2, 


1442.2] 


[1388.5 


1460.6] 


-1.09 





610±0 


031 


V*BNPav 


2283.5 


[2263 


4. 


2303.8] 


[2243.4 


2324.2] 


0.89 





511±0 


026 


V*BNVul 


680.4 


[674 


6, 


686.4] 


[668.7 


692.4] 


-1.99 





819±0 


039 


V*BPPav 


2248.3 


[2231 


5, 


2265.3] 


[2214.8 


2282.3] 


0.00 





587±0 


025 


V*CGLib 


1091.5 


[1073 


7, 


1109.6] 


[1056.2 


1128.0] 


0.12 





638±0 


059 


V*CIAnd 


1899.9 


[1885 


2, 


1914.8] 


[1870.6 


1929.7] 


0.65 





683±0 


021 


V*CNLyr 


1065.1 


[1055 


0, 


1075.3] 


[1044.9 


1085.6] 


0.23 





764±0 


032 


V*DDHya 


2011.6 


[1996 


2, 


2027.2] 


[1980.9 


2042.9] 


0.41 





651±0 


023 


V*FWLup 


409.7 


[406 


7, 


412.7] 


[403.7 


415.7] 


1.34 





696±0 


019 


V*HHPup 


1054.3 


[1043 


3, 


1065.5] 


[1032.5 


1076.7] 


0.99 





679±0 


030 


V*HKPup 


1294.6 


[1277 


5, 


1311.9] 


[1260.6 


1329.4] 


3.66 





192±0 


036 


V*IOLyr 


1601.8 


[1588 


4, 


1615.4] 


[1575.1 


1629.0] 


0.38 





618±0 


024 


V*MSAra 


1621.3 


[1608 


9, 


1633.8] 


[1596.5 


1646.4] 


-0.41 





639±0 


021 


V*MTTel 


429.7 


[423 


1, 


436.4] 


[416.6 


443.2] 


-1.35 





678±0 


035 


V*RRGem 


1242.5 


[1229 


8, 


1255.4] 


[1217.2 


1268.3] 


1.26 





687±0 


025 


V*RRLyr 


252.9 


[250 


9, 


254.9] 


[248.9 


256.9] 


-1.16 





636it0 


018 


V*RSBoo 


840.2 


[830 


7, 


849.8] 


[821.3 


859.6] 


0.93 





718±0 


026 


V*RVCct 


1070.8 


[1060 


1, 


1081.6] 


[1049.6 


1092.5] 


-0.98 





680±0 


023 


V*RVCrB 


1313.1 


[1294 


3, 


1332.2] 


[1275.7 


1351.5] 


-1.18 





692±0 


032 


V*RVOct 


984.3 


[976 


3, 


992.4] 


[968.4, 


1000.5] 


0.50 





469±0 


040 


V*RWCnc 


1704.9 


[1691 


3. 


1718.6] 


[1677.8 


1732.4] 


-0.67 





629±0 


024 


V*RWDra 


1579.6 


[1566 


5, 


1592.9] 


[1553.5 


1606.2] 


-1.14 





716±0 


020 


V*RWTrA 


1122.5 


[1109 


5, 


1135.6] 


[1096.7 


1148.8] 


0.14 





879±0 


027 


V*RXCol 


2272.3 


[2252 


3, 


2292.4] 


[2232.5 


2312.7] 


-0.78 





637±0 


029 


V*RXEri 


590.1 


[585 


1, 


595.1] 


[580.2 


600.1] 


-0.06 





630±0 


020 


V*RYCol 


1086.1 


[1078 


0, 


1094.3] 


[1070.0 


1102.5] 


0.53 





651±0 


018 


V*RYOct 


1761.4 


[1747 


2, 


1775.7] 


[1733.0 


1790.2] 


0.48 





446±0 


022 


V*RZCep 


381.3 


[375 


2, 


387.6] 


[369.2 


393.9] 


-0.99 





650±0 


038 


V*RZCet 


1604.0 


[1591 


8, 


1616.3] 


[1579.6 


1628.8] 


-0.63 





695±0 


020 


V*SAra 


941.0 


[933 


3, 


948.8] 


[925.6 


956.6] 


1.60 





555±0 


020 


V*SSOct 


1144.8 


[1133 


9, 


1155.8] 


[1123.1 


1167.0] 


-0.07 





569±0 


030 


V*STBoo 


1188.8 


[1177 


2, 


1200.5] 


[1165.7 


1212.3] 


-0.38 





572±0 


022 


V*STVir 


1468.1 


[1453 


9, 


1482.4] 


[1439.8 


1496.9] 


1.59 





564±0 


024 


V*SUDra 


696.2 


[688 


8, 


703.7] 


[681.4 


711.3] 


0.01 





513±0 


025 
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Table 2 — Continued 



]\T Q TVl O 

IN Ojillt: 






fj 


d + la] 


[d — 2 


cr, d + 2o-] 






Mv 


C 




Inr] 




[pc] 




IP 




fNn nf al 




[mag] 


V*SVEri 


702.4 


[693 


u 


Vl 1 ol 
( ii.Zj 


[685 





720.2] 


0.91 





A on J_n 


UoU 


V*SXFor 


1197.5 


[1186 


O 


izUo.Oj 


[1175 


9 


1219.5] 


-1.08 





DooitU 




V*SZGem 


1622.1 


[1609 


/ 


iDo4.Dj 


[1597 


5 


1647.1] 


0.27 





oooitu 


UzU 


V*TTCnc 


1227.2 


[1217 


q 
O 


iZo ( .oj 


[1207 


5 


1247.3] 


-0.87 





Dm itU 




V*TTLyn 


667.7 


[661 


O 


1 O. i \ 


[656 





679.7] 


-0.92 





DoDitU 




V*TVCrB 


1764.9 


[1749 


/ 


1 7sn ol 
1 ( ttU.zj 


[1734 


7 


1795.6] 


-1.00 







UZz 


V*TWHer 


1246.4 


[1233 


Q 
O 


1 oc;q 1 1 
izoy . ij 


[1221 


4 


1271.9] 


0.79 





DOOltU 




V*TWLyn 


1677.6 


[1664 


7 


loyu.Dj 


[1651 


9 


1703.7] 


0.34 





I ozitu 




V*TYAps 


1422.3 


[1411 


Q 
O 


1 /I QO ol 

i4oz.yj 


[1401 


3 


1443.6] 


0.35 





DD4itU 




V*TZAur 


1605.2 


[1588 


2 


1622.4] 


[1571 


3 


1639.8] 


0.44 





689±0 


027 


V*ULep 


976.8 


[968 


6 


985.1] 


[960 


4 


993.5] 


0.02 





516±0 


020 


V*UPic 


1357.3 


[1346 


1 


1368.6] 


[1334 


9 


1380.0] 


0.45 





703±0 


021 


V*UVOct 


535.8 


[531 


8 


539.8] 


[527 


8 


543.9] 


-0.56 





598±0 


017 


V*UYBoo 


1154.2 


[1141 


9 


1166.6] 


[1129 


8 


1179.1] 


-1.53 





542±0 


024 


V*UYCam 


1442.5 


[1413 


8 


1471.7] 


[1385 


6 


1501.6] 


-0.32 





665±0 


044 


V*V413CrA 


852.3 


[844 


9 


859.8] 


[837 


5 


867.3] 


-0.33 





680±0 


022 


V*V440Sgr 


770.5 


[764 


6 


776.5] 


[758 


8 


782.5] 


-0.09 





617±0 


020 


V*V4450ph 


773.5 


[765 


5 


781.5] 


[757 


7 


789.6] 


1.72 





640±0 


037 


V*V4550ph 


1816.7 


[1801 


8 


1831.8] 


[1787 





1847.0] 


0.78 





581±0 


026 


V*V499Cen 


1133.1 


[1124 


2 


1142.1] 


[1115 


4 


1151.1] 


0.06 





591±0 


021 


V*V675Sgr 


803.2 


[795 


1 


811.5] 


[787 





819.8] 


0.39 





353±0 


026 


V*VInd 


698.9 


[693 


6 


704.2] 


[688 


3 


709.6] 


-0.03 





587±0 


018 


V*VXHer 


970.8 


[963 





978.7] 


[955 


2 


986.7] 


-0.25 





596±0 


020 


V*VYLib 


1256.2 


[1246 


2 


1266.2] 


[1236 


3 


1276.3] 


-0.17 





642±0 


027 


V*VYSer 


768.6 


[759 


1 


778.3] 


[749 


6 


788.1] 


-0.41 





569±0 


028 


V*VZHer 


1442.2 


[1430 





1454.4] 


[1417 


9 


1466.8] 


0.77 





594±0 


020 


V*WYPav 


1641.2 


[1627 





1655.7] 


[1612 


8 


1670.2] 


-0.58 





778±0 


022 


V*XAri 


500.5 


[495 


3 


505.8] 


[490 


2 


511.1] 


-0.52 





438±0 


032 


V*XXAnd 


1030.7 


[1017 


5 


1044.1] 


[1004 


4 


1057.6] 


0.05 





476±0 


029 


V*XXPup 


1225.1 


[1215 


9 


1234.4] 


[1206 


7 


1243.8] 


0.16 





601±0 


019 


V*XZAps 


1880.7 


[1864 


6 


1897.1] 


[1848 


5 


1913.5] 


0.91 





565±0 


031 



''Best distance posteriors from the analysis described in §4. 

''The number of a discrepancy between the prior and posterior mean, defined as A(/iprior — 

-1 Mprior ~ Apost 

/ipostj — r-2 —2 

V'"f'prior+%po=t 

'^Absolute V-band magnitudes calculated with posterior from the PLR analysis (§1) and 
from converted, extinction corrected Hipparcos mean-magnitude, (§3; Table 1). 



